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We propose a superconducting qubit design, based on a tunable RF-SQUID and nanowire kinetic 
inductors, which has a dramatically reduced transverse electromagnetic coupling to its environment, 
so that its excited state should be metastable. If electromagnetic interactions are in fact responsible 
for the current excited-state decay rates of superconducting qubits, this design should result in 
a qubit lifetime orders of magnitude longer than currently possible. Furthermore, since accurate 
manipulation and readout of superconducting qubits is currently limited by spontaneous decay, 
much higher fidelities may be realizable with this design. 

PACS numbers: 



One of the distinguishing features of Josephson- 
junction (JJ)-based qubits is their strong couphng to 
electromagnetic (EM) fields, which permits fast gate op- 
erations (~ 10-lOOns). However, it may also be re- 
sponsible for their short excited-state lifetimes (< 4/is 
0)0)11) 0,111); that is, assuming the decay process is elec- 
tromagnetic, its rate depends on the same matrix element 
which governs intentional qubit manipulations by exter- 
nal fields. Unfortunately, understanding and controlling 
spontaneous decay of these circuits has so far proved diffi- 
cult, because it also depends on their EM environment at 
GHz frequencies. This environment is affected not only 
by packaging and control lines, but also by microscopic 
degrees of freedom in the substrate, surface oxides, and 
JJ barrier dielectrics. In fact, low-frequency noise due 
to microscopic fluctuators is already known to produce 
"dephasing" of qubits P, 0, H, 0, [1] . Although little is yet 
certain about the properties of these degrees of freedom, 
work is ongoing to study them and to reduce their 
number through improved materials and fabrication 

In this Letter we discuss a different approach, seeking 
a qubit which is insensitive to high-frequency EM fluctu- 
ations by design. This is a departure from the research 
area known as circuit QED [8|, in which strong trans- 
verse coupling to EM fields is both a prerequisite and a 
figure of merit. We will show that a qubit design based 
on weak transverse EM coupling could yield significantly 
longer excited-state lifetimes, while still allowing practi- 
cal, scalable computation. 

The decay rate of an excited state |e) to lower-lying 
state \g) is typically given by Fermi's golden rule: F = 
l/Ti = (2TT/h)\mi\'^p{hweg), where rm = {e\Hi\g), Hi is 
a Hamiltonian describing the coupling between the qubit 
and a continuum (e.g., the excited states of an ensemble 
of two- level systems - TLSs), and pihweg) is the density of 
states in that continuum at the energy iTiHeg = E,, — Eg. 
The rui can be nonzero for a JJ-based qubit when the 
fiux through a loop, the induced charge across a JJ, or a 
JJ critical current depends on the state of one or more 
TLSs. To minimize the resulting decay rate, we must 
reduce p ov mi. Our focus here will be on the latter. 

A good choice for qubit energy levels which are weakly 
coupled to each other by EM fields are the flux states 



of an RF SQUID [Fig. [^a)] [S] at large Ej/Ec, where 
Ej — <i>o-fc/27r and Ec = jlCtot are the Josephson 
and charging energies {Jc is the J J critical current, $o = 
hjle, and Ctot is the total capacitance across the JJ). 
When ^ 'fo/2, two quantum states, in which either 
zero or one fluxon is contained in the loop, become nearly 
degenerate, and separated by a potential barrier [Fig. 
mb)]. The Hamiltonian for the RF SQUID is ;1Q]: 

H ^ AEc{n~ Uef - Ejcosij) + Eli'^ /2 (1) 

where (j) is the phase across the J J, h = —id/dcj) is 
operator corresponding to the number of Cooper pairs 
that have tunneled through the J J, 7 = (f) + 2nf is the 
phase across the inductor, / — ^rf/^q^ and El = 
(<I'o/27r)^/i. The quantity is a fluctuating offset 
charge across Ctot induced by capacitances to the en- 
vironment or by tunneling of quasiparticles through the 
junction (at DC ne=0 due to the inductive shunt). 

We diagonalize H on a. lattice of 7 points to obtain 
wavefunctions ^k{l) = (7!^) [Fig. H^c), (e), and (g)], 
which are then used to evaulate transition matrix el- 
ements {k\Hi\k') [ll| for flux, charge, and /c-coupled 
TLSs, with: Hf = 2'KSfEjsm{'y + 27r/), H„ = SSnEcfi, 
and Hi = SIcE j cos{j + 2TTf), respectively, and 6f, 
Sn, and Sic are the (small) amplitudes of TLS-state- 
depcndent changes in /, rig, and Ic- These amplitudes 
will be different for each TLS, so it is conceptually use- 
ful to recast the golden rule in terms of an average noise 
power spectral density Si jTlJ thus: F^ = \di\'^ Si{u!eg)/h'^ , 
where di = (e|Xi|g) are analogous to a transition dipole 
for each fluctuation, and: Xf = 27r£^jsin(7 + 27r/), 
Xn = 8Ecn, Xi^ = Ej 005(7 -|-27r/) with units of energy 
per $0) electron pair, and current. 

Since the operators Xi are local in 7, a simple way to 
reduce all of the di at once is to reduce the overlap of 
the probability distributions IV'glT)!^ and |^e(7)|^. This 
overlap results from tunneling through the barrier [Fig. 
[IJb),(c)], so to minimize it we detune the left and right 
wells from each other (/ ^ 0.5) to avoid resonant tunnel- 
ing, and increase the barrier height by increasing Ej / Ec 
and/or Ej/El [Fig. [Ii:d),(e)]. 
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FIG. 1: (color online) Fluxon tunneling for the RF SQUID flux qubit. (a) schematic, (b) potential, and (c) qubit level 
wavefunctions \g) (solid) and |e) (dashed) at / = 0.515 for Ej, Ec, El — hx 120, 6, 60 GHz. (d) potential, and (e) wavefunctions, 
for Ej,Ec,El = h X 180,4,60; the potential barrier between wells is higher so the tunneling is weak, (f) potential, and (g) 
wavefunctions for Ej,Ec,El — hx 120,6,0.375. Dotted lines in (c), (e) and (g) show the next excited states; In (c) and (e) 
these are "vibrational" excitations, while in (g) they are the ground states of adjacent wells (corresponding to -1 or 2 fluxons 
in the SQUID loop; fluxon numbers for each well are indicated in (f)). 



Unfortunately, when / 7^ 0.5, dujpg/df ^ 0, and 
nonzero, low- frequency 5 J produce dephasing ilj. This 
sensitivity can be reduced by increasing L, since hiiOeg ~ 

^(/— 0.5), for El ^ Ej [l^. To realize large L, increas- 
ing the loop size is not attractive, both because it would 
need to be of millimeter scale, and because its large ca- 
pacitance would limit Eq- Instead, we propose using the 
kinetic inductance of a long, meandered nanowire pat- 
terned from thin (^5 nm thick) NbN, which can have 
sheet inductance as large as ~100 pH and Ic ^ 20/iA 
[l^. [l5|. A 10 /xm-square meander of 100 nm-wide wire 
gives L ~500 nH [l^, and EM simulation shows a shunt 
capacitance of only ~0.4 fF, significantly smaller than 
that of the JJs we consider below (~3.2 fF). 

Figure [2l[a) shows the resulting \di\ for our proposed 
qubit, as a function of Ej/Eq- Also shown, by horizon- 
tal dashed lines, are the \di\ for transmon, quantronium, 
flux, and phase qubits Jjj. Based on these results, and 
by extracting bounds on the Si from Ti values observed in 
Refs. we can estimate Ti for our qubit. Not 

surprisingly, no single set of 5*^, in conjunction with the 
calculated di, can accurately explain all of the observa- 
tions, since the noise levels are likely somewhat different 
in each experiment; however, for the present purpose, we 
take: Si{5GEz) < 1.4 x 10-iVA^Hz-\ from Ti = 650ns 
for the phase qubit of Ref. Q; 5'„(5.7GHz) < 1.6 x lO^^^ 
Hz^^ from Ti = 1.7 fis for the transmon of Ref. 0]; and 
S'/(5.5GHz) < 1.3 X 10-20 Hz"! from Ti = 2^s from the 
flux qubit of Ref. [l|. Panel (b) shows the resulting es- 
timate of Ti for our qubit (dominated by charge noise). 
For Ej/Ec ~ 3 (as in Ref. [H), Ti ~ 3fis, similar to 
what was observed; however, at Ej/Eq — 20, we find 
Ti - 950 milliseconds 

The reduced transverse coupling that we achieve 
through increasing Ej / Eq also means we must drive the 
qubit with larger fields to manipulate it [l8|, or mea- 
sure it dispersively [H, [l^ . At some point the required 
driving becomes strong enough that spurious effects oc- 
cur, such as off-resonant excitation of strong transitions 
to short-lived excited states (followed by decay) , or large 
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FIG. 2: (color online) Transverse coupling of the metastable 
RF SQUID qubit vs. Ej/Ec, for Ec,El ^ hx6, 0.375 GHz, 
and / = 0.57. Panel (a) shows Ke|X„|g)| (dashed), \{e\Xf\g)\ 
(dash-dot), and |(e|X/|(;)| (solid), respectively. Horizontal 
lines show equivalent |(e|Xi|(;)| for the transmon quantro- 
nium [^, flux P, 0], and phase [3| qubits. (b) shows the pre- 
dicted excited state lifetime for the metastable RF SQUID 
qubit. (c) schematic for the proposed qubit, with Ej/Ec 
tunable through the flux $dc. The gradiometric design de- 
couples this flux from e = /p($i — $2) [23 |. 

nonlinearities in the qubit response. Furthermore, ini- 
tializing the qubit state will take longer as the lifetime 
is increased [2l| . It will therefore be useful to be able to 
adjust Ej/Ec in real time using a tunable RF SQUID 
[Fig. [He)] (analogous to the tunable flux qubit [13, H^). 
The single JJ is replaced by a DC SQUID, and the RF 
SQUID loop is replaced with a gradiometric design where 
fjiF = ($1 — ^2)/^o- In this configuration, Ej in eq. [T] 
is replaced with: 

Ejifoc) = 2E JO coa[7: f DC -1 DC m (2) 
« 2E'joCoa[7:fDc] (3) 

where Ejq is the Josephson energy of each J J, foe = 
^Dc/^o, and ^dc is the phase across Ljjc, the self- 
inductance of the DC SQUID loop. To obtain eq. [31 
we note that for L^c ^ L,Lj = $o/27r/c, the zero- 
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FIG. 3: (color online) Transverse manipulation of the 
metastable RF SQUID qubit. Integration of the Master equa- 
tion (including decay) for the first 10 levels of the RF SQUID, 
undergoing a 7r-pulse starting from \g). The solid line is for 
the modulation of <1>_dc in the left inset; the dashed line is for 
a sinusoidal modulation of Vrf with amplitude We = 0.24, at 
E,j = /ix42 GHz. The dotted line is the equivalent result for 
a flux qubit 



point fluctuations oi^oc can be adiabatically eliminated, 
yielding to leading order only a small renormalization of 
EjQ [24, ^25] {ior Ldc < 50 pH, and the parameters under 
consideration here, a fraction of a percent). 

The qubit can be manipulated (or measured disper- 
sively) with Vrf, $dc, or <^rf = $1 - $2 [Fig. HJc)]. 
We consider the first two here. In order to describe large- 
amplitude driving, and to incorporate spontaneous decay 
between instantaneous energy eigenstates |m(t))', we use 
a time-dependent transformation to the instantaneous 
energy eigenbasis, yielding the Hamiltonian: Had — 
RHW — ihR^W where H is given by eq. [T] and R 
is defined by: R\ip) = \tp') (the prime indicates the ba- 
sis |r7i(i))'). The first term is diagonal, containing the 
time-dependent eigenenergies, and the second term yields 
transitions between levels. We integrate a master equa- 
tion based on Had, truncated to the 10 lowest-lying in- 
stantaneous eigenstates (up to >100 GHz above \g)). To 
this we add a spontaneous decay rate 7mn(0 from each 
level \m') to each other level |n'). To generate the 7mn(i), 
we use Fermi's golden rule, and assume an Ohmic noise 
spectrum Si{Lu) oc huj/{l — p-^^/^sT-^ (w > denotes 
downward transitions, cj < upward), with the overall 
amplitude for each type of noise discussed above. The 
time dependence of the 7m„(t) comes from the \Xi{t)\'^^^. 

As a test case, we consider a 7r-pulse, where the qubit 
starts in \g), for which an indication of gate fidelity is 
how much population we can put in |e), as shown in Fig. 
H We take Ec, El ^ hx6, 0.375 GHz (L = 430 nH), and 
Jrf = 0.57 (cJeg = 27rx 1.034 GHz For modulation 

of $ DC (solid line) , we use the pulse shown in the left in- 
set to Fig. [3l which starts and ends at Ej = hx 200GHz 
(with E'jg — hx280 GHz). For modulation of Vrf, we 
take a fixed Ej = hx42 GHz. The simulation yields 
1 - Pe = 1.1 X IQ-s and 1 - Pe = 2.5 x IQ-^ for ^rf 



and Vrf modulation, respectively. The former is limited 
almost completely by decay of |e) during the brief ex- 
cursions to smaller Ej/Eq where 710 (i) is larger. This 
also explains the shape of the time evolution: the drive 
becomes effectively faster when Ej/Ec is smaller, pro- 
ducing the upward "steps" . Spurious excitation to adja- 
cent fiuxon states (-1,2 in Fig. [IJf)) and higher vibra- 
tional states is at the O(10~^) level. Driving with Vrf is 
limited by off-resonant excitation of the first vibrational 
levels (at ~40 GHz) followed by decay. This process is 
suppressed for $ dc modulation since the perturbation is 
nearly even about the potential well center. For compar- 
ison is shown the same simulation for a flux qubit [If, 0] , 
which has 1 — Pe = 2 x 10~^, due to decay from |e). 

This simulation does not include 1/f flux noi se |l| . To 
estimate its effect, we use the results of Ref. [231, ^^i*^ 
the fact that for L = 430 nH, dcjeg/df = 14.3 MHz/m$o 
(^100 times smaller than a typical flux qubit far from 
/ — 0.5). For the noise amplitude measured in Ref. [l|, 
we calculate the average error in the qubit relative phase 
over the 8 ns 7r-pulse to be ^ 4.5 mrad, which for the 
maximally sensitive (l^) -f- |e))/-\/2 state gives an error 
probability of only ~ 2.0 x 10"^ 

By eliminating the transverse coupling between qubit 
levels induced by external fields, we have also eliminated 
the usual mechanism for coupling qubits to each other 
[29| . Instead, we must use a longitudinal coupling, sim- 
ilar to Ref. [30|, which makes use of the nonzero flux 
tunability [3l|. A schematic of our proposed circuit is 
shown in Fig. (Ha). Two RF SQUID qubits are cou- 
pled by mutual inductances M to a third coupler qubit 
with large persistent current Ip , biased at its degeneracy 
point {Jrf = 0.5). The approximate Hamiltonian is: 

1,2,C 1,2 
i j 

Here, eigenstates of ct^ are well-defined fiux (persistent 
current) states, and Jc = MI^ dei^2/d'^ with dti^^/d'^ « 
Att'^El/^o m, Jo = Molplp. The Sj are residual trans- 
verse flux coupling of the data qubits due to nonzero df. 
We take: = 5.2/zA and tc = hx 5GHz M = 5pH, 
Mo = O.lpH, L = 430nH, and I^-^ « $o/2P = 2.4 nA, to 
obtain Jc = /ixl88 MHz, Jq ^ h x 0.87 kHz; this gives 
a conditional frequency shift h6i, ~ 2J^/tc — Jq = 14.1 
MHz ^] and a conditional- 7r/2 gate in « 18 ns. If we use 
spin-echoes [13, HO] during the gate, the residual phase 
drift due to 1 /f flux noise (not canceled by the echo) dur- 
ing this time is 3.6 mrad, producing a maximal error (in 
addition to that from the tt pulse) of 3.0 x 10"'' [28| . 

A very small transverse coupling to the data qubits 
{Sj <C 1) also means that their excited states will un- 
dergo negligible mixing with the excited state of the cou- 
pler (which will likely be short-lived). Figure Ub) shows 
the decay rates that result. These are proportional to: 
|(zj|(Tp|/cZ)P {i,j,k,l G {g,e}), where are the com- 
putational states (the lowest four eigenstates of eq. SI 
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FIG. 4: (color online) Switchable coupling between 
metastable RF SQUID qubits. (a) schematic; two data 
qubits (RF SQUIDs with meandered kinetic inductors) are 
coupled through a mutual inductance M to a coupler qubit 
(RF SQUID with large Jp, and ec = 0) and parasitically to 
each other (Mo), (b) shows the calculated decay rates for 
the four computational levels, relative to the decay rate of 
the coupler in isolation, for the parameters in the text, and 
with e\,€2 — h X 1.0, l.lGHz (chosen to be different only for 
clarity). The inset shows the resonances that occur due to 
nonzero Si. 



which in the Jc limit correspond to the coupler in 
its ground state [13]). The pronounced peaks (and dips) 
occur when the coupler is nearly resonant with one of 
the data qubits; in these regions, the nonzero Sj produce 
two entangled states of a data qubit and the coupler, 
with one state coupling maximally to fluctuations and 



the other minimally. When both qubits are detuned far 
from the coupler, the data qubit decay rate is sufficiently 
suppressed that even coupler qubit lifetimes at the ns 
scale would have little effect. 

In summary, we have described a qubit design with 
weak transverse coupling to EM fields. This qubit should 
be significantly less sensitive to microscopic EM degrees 
of freedom arising from fabrication imperfections, and 
may permit very long Ti times with good fabrication 
yield. Although these qubits are still weakly sensitive 
to low-frequency flux drifts, this is in principle a problem 
with any manufactured qubit; e.g., even a flux qubit at its 
flux-insensitive point is sensitive to drifts of Ic at a level 
that may soon become a coherence limit 33|, and even 
lithographically defined superconducting resonators have 
resonance- frequency noise [s^ . Thus, man-made qubits 
may inevitably require modified encoding/computation 
schemes which are resistant to the inevitable drift be- 
tween each qubit 's relative phase and the absolute phase 
reference required for sustained computation. 

We acknowledge helpful discussions with William 
Oliver, Jeremy Sage, and Jens Koch. 

This work is sponsored by the United States Air Force 
under Contract #FA8721-05-C-0002. Opinions, inter- 
pretations, recommendations and conclusions are those 
of the authors and are not necessarily endorsed by the 
United States Government. 



[1] F. Yoshihara, K. Harrabi, A. O. Niskanen, Y. Nakamura, 
and J. S. Tsai, Phys. Rev. Lett. 97, 167001 (2006). 

[2] P. Bertet, et al, Phys. Rev. Lett. 95, 257002 (2005). 

[3] J. A. Schreier, et al, Phys. Rev. B 77, 180502(R) (2008). 

[4] M. Hofheinz, et al, Nature 459 5463 (2009). 

[5] 1. Siddiqi, et al, Phys. Rev. B 73, 054510 (2006). 

[6] Z. Kim, et al, Phys. Rev. B 78, 144506 (2008); S. Sendel- 
bach, et al, Phys. Rev. Lett 100, 227006 (2008); R. W. 
Simmonds, et al, Phys. Rev. Lett. 93, 077003 (2004); A. 
Lupascu, P. Bertet, E.F.C. Driessen, C.J. P.M. Harmans, 
and J.E. Mooij, arXiv 0810.0590 (2008). 

[7] J. M. Martinis, et al, Phys. Rev. Lett. 95, 210503 (2005); 
J.S. Kline, H. Wang, S. Oh, J.M. Martinis and D.P. Pap- 
pas, Supercond. Sci. Tech. 22 015004 (2009); S. Oh, et 
al, Phys. Rev. B 74, 100502(R) (2006). 

[8] R.J. Schoelkopf and S.M. Girvin, Nature 451, 664 (2008), 
and references therein. 

[9] J.R. Friedman, V. Patel, W. Chen, S.K. Tolpygo, and 
J.E. Lukens, Nature 406, 43 (2000). 
[10] J. Koch, V. Manucharyan, M.H. Devoret, and L.I. Glaz- 

man, arXiv:0902.2980 (2009). 
[11] J. Q. You, X. Hu, S. Ashhab, and F. Nori, Phys. Rev. B 

75, 140515 (2007). 
[12] An alternative is a series array of JJs, as demonstrated 
in Ref. [l^l . However, since each of these JJs is a tunnel- 
ing path for the fluxon, this approach cannot be used 
to strongly suppress di in the manner discussed here. 
Although fluxons can in principle tunnel through NbN 



nanowires, due to their granular microstructure or quan- 
tum phase slips, we estimate these processes occur at 
rates < 0(10 s"^). 

[13] V. Manucharyan, J. Koch, L.I. Glazman, and M.H. De- 
voret, arXiv:0906.0831 (2009) 

[14] Note that Ref. ^Qj was in the opposite limit Ej <C El- 

[15] A.J. Kerman, et al, Appl. Phys. Lett. 88, 111116 (2006). 

[16] For the transmon and quantronium, the Ti is only sen- 
sitive to charge noise, and for the flux qubit at / = 0.5, 
charge and flux noise. For all three, we diagonalize the 
Hamiltonian in the charge basis. 

[17] One possible non-EM decay mechanism is acoustic cou- 
pling of an oscillating persistent current to phonons [E.M. 
Chudnovsky and A.B. Kuklov, Phys. Rev. B 67, 064515 
(2003)]; however, the oscillation frequency is limited by 
the tunelling rate through the E,j barrier. 

[18] Raman transitions through an excited state have also 
been proposed [M.H. Amin, A. Yu. Smirnov, and A. 
Maassen van den Brink, Phys. Rev. B 67, 100508(R) 
(2003)1; however, the decay rates of relevant higher states 
may be too large to do this effectively. 

[19] Weak coupling (large E,j/Ec) may also be interesting 
for dispersive readout. In contrast to the strong-coupling 
limit, in which dephasing due to photon number splitting 
is important for weak coupling large fields could 

be used, with a corresponding reduction in (fractional) 
photon-number uncertainty. 

[20] J. Gambetta, et al, 042318 (2006). 



5 



[21] One might use active cooling, as in [S.O. Valenzuela, et 
al, Science 314, 1589 (2006)] rather than waiting for 
the qubit to decay; however, this wiU also become more 
difficult at larger Ej/Ec- 

[22] T. P. Orlando, et al, Phys. Rev. B 60, 15398 (1999). 

[23] F. G. Paauw, A. Fedorov, C. J. Harmans, and J. E. 
Mooij, Phys. Rev. Lett. 102, 090501 (2009). 

[24] A. Maassen van den Brink, Phys. Rev. B 71, 064503 
(2005). 

[25] The 7_DC oscillation frequency is ~ C(THz), so yoc 
can be replaced by the ground state probability distri- 
bution: P(7Dc) = (2/7r7jJ,)exp[-2(7i3c - {7dc>)^/7^]i 
with 7„ = 2n{e^LDc/^oCjy^\ (jdc) = Ejo/EE^[1 - 
cos(7r/i3c)l < 7uj- Averaging eq. [2] over P{'yDc) gives: 
Ejifoc) ~ Ejoexp[~'yl/128]. 

[26] A relatively small uieg is used to avoid spurious excitation 
of the adjacent fluxon states by a short pulse (at L = 
430nH these states are only ~ 15 GHz away). 

[27] G. Ithier, et al, Phys. Rev. B 72, 134519 (2005). 

[28] These phase drifts continue to accrue with time, even 
when spin-echoes are used (e.g., eq. 22 of Ref. [13]). For 



the present qubit, this is (13/^s) ^ for the flux noise mea- 
sured in Ref Q. 

[29] A.O. Niskanen, et al, Science 316 723 (2007); T. Hime, 
et al, Science 314 1427 (2006); S.H.W. van der Ploeg, 
et al, Phys. Rev. Lett. 98, 057004 (2007); L. DiCarlo, et 
al, Nature 460 240 (2009); S. Ashhab, et al, Phys. Rev. 
B 77, 014510 (2008). 

[30] A.J. Kerman and W.D. Oliver, Phys. Rev. Lett. 101, 
070501 (2008). 

[31] A related proposal was detailed in [Y.-D. Wang, A. 

Kemp, and K. Semba, Phys. Rev. B 79 024502 (2009)]. 

In this scheme the coupling is switched off by tuning the 

qubits to denegeracy, so it would not be applicable here. 
[32] This 7^ results from: Ej = 5000 GHz, Ec = 0.7 GHz, 

L = 35 pH. 

[33] D.J. Van Harlingen, et al., Phys. Rev B 70, 064517 
(2004). 

[34] R. Barends, et al., IEEE Trans. Appl. Sup. 19, 936 
(2009). 



